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Abstract 

Recently, the HESS telescope has observed a gamma-ray spectrum from the galactic center extending 
to at least ~10 TeV. Although it has been suggested that this could be the product of annihilating 
dark matter particles, the candidates most frequently discussed (such as neutralinos) are far too light 
to account for this flux. In this letter, we consider stable particles from the messenger sector of gauge 
mediated supersymmetry breaking models as an alternative dark matter candidate. We find that a 
20 to 30 TeV messenger state can provide a thermal relic density consistent with the measured dark 
matter density of the universe and can indeed generate the spectrum observed by HESS. 
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I. INTRODUCTION 



Recently, the High-Energy Spectroscopic System (HESS) collaboration has reported the 
detection of gamma-rays from the region of the galactic center with energies ranging between 
approximately 250 GeV and 10 TeV [l[. Two other Atmospheric Cerenkov Telescopes (ACTs), 
Cangaroo-II [H and Whipple [j| have also reported positive detections of very high-energy 
gamma-rays from the region, although with considerably lower significance. 1 

Currently, the origin of this very high-energy emission is unknown. Efforts have been made 
to explain the data with astrophysical processes involving the supermassive black hole known 
to exist at the galactic center |4|, |5J]. Although these are certainly possible explanations, the 
source of the observed emission is not yet clear. 

Gamma-ray observations of the galactic center have long been studied as one of the most 
promising methods to search for particle dark matter {gj. The most popular dark matter 
candidates 0], such as the lightest neutralino in _R-parity conserving supersymmetric models, 
are generally much lighter than the highest energy gamma-rays seen by HESS, and therefore 
cannot possibly generate the observed spectrum. It has been shown that a dark matter particle 
annihilating through typical channels (qq, W + W~ , Z°Z°) would require a mass between 12 and 
29 TeV to generate this spectrum 0. In addition to the unacceptable levels of fine tuning 
which would be required, a neutralino in this mass range would generate far more dark matter 
during thermal freeze-out than is observed. For these reasons, we do not consider a very heavy 
neutralino to be a likely source of the flux observed by HESS. 

Models of Gauge Mediated Supersymmetry Breaking (GMSB) provide an interesting alter- 
native to neutralino (or other lightest supersymmetric particle) dark matter. Typically in these 
models, supersymmetry breaking originates in a strongly interacting sector (the supersymme- 
try breaking sector) and is communicated to the superpartners of the Standard Model (SM) 
states by SM gauge interactions with particles within a messenger sector It is natural for 
the lightest messenger particle to be stable, electrically neutral and have a mass on the order 
of a few times 10 TeV. Such a particle, as we will show, can have an annihilation cross section 
which generates a thermal relic density consistent with the measured dark matter density of the 
universe. In this letter, we calculate the gamma-ray spectrum generated in the annihilations 
of GMSB messenger dark matter and show that this particle can be a viable source of the very 
high-energy gamma-rays observed by HESS. 

II. MESSENGER DARK MATTER IN GAUGE MEDIATED SUPERSYMMETRY 
BREAKING MODELS 

GMSB models typically contain a messenger sector, charged under the SM gauge group, 
which communicates supersymmetry breaking to the superpartners of SM states by SM gauge 
interactions. Such gauge communication of supersymmetry breaking has the great advantage 
that unobserved flavour-changing-neutral-current, CP-violating, and rare decay processes aris- 
ing from the soft-supersymmetry-breaking terms are automatically absent, in contrast to most 
other methods of communication. Since supersymmetry breaking is communicated by renormal- 



The spectrum reported by the Cangaroo collaboration differs substantially from the HESS results. The results 
of the Whipple telescope, which are limited to an integrated flux, are only marginally consistent with HESS. 
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izable gauge interactions, the fundamental scale of supersymmetry breaking in GMSB models is 
much lower than for supergravity communication, possibly as low as 50 to 100 TeV depending 
on the interactions between the superymmetry breaking and messenger sectors. As a conse- 
quence, the Lightest Supersymmetric Particle (LSP) is the gravitino, with mass, m 3 / 2 ~ 1 — 10 
eV. Such a mass is far too small to contribute substantially to the dark matter density of the 
universe. Thus, in this class of models, it is necessary to look beyond the LSP for the con- 
stituents of dark matter. Fortunately in the simplest class of GMSB models, the messenger 
sector naturally contains a stable state with mass of a few ten's of TeV. 2 

Messenger particles in GMSB models have been studied as potential dark matter candidates 
in the past. In Refs. 11, minimal supersymmetric models with gauge mediated supersym- 



metry breaking were studied. In these models, the messenger sector consists only of particles 
carrying SU(3)xSU(2)xU(l) charge, the lightest of which will be stable if the supersymmetry 
breaking sector contains only singlets charged under this group. In substantial portions of the 
parameter space, the lightest messenger state is a massive complex scalar with the SM gauge 
quantum numbers of a left handed neutrino. 

This simple model ultimately fails to provide a viable dark matter candidate, however. With 
only gauge interactions to contribute to its annihilation rate in the early universe, it has been 
shown that such a particle must be lighter than a few TeV to not overdose the universe (ll. 13 . 



In addition to being too light to produce the very high-energy gamma-rays observed by HESS, 
such a dark matter candidate is excluded by direct detection searches. 

A very simple way to extend the Minimal Supersymmetric Standard Model (MSSM) is 
to introduce an additional Higgs singlet. This model, the Next-to-Minimal Supersymmetric 
Standard Model (NMSSM), has several advantages compared to the MSSM. Most notably, it 
provides a natural mechanism for the /^-parameter to obtain a value near the electroweak scale 
(which is required for acceptable electroweak symmetry breaking), rather than Mqvt or Mp\. 
It also ameliorates the mild fine tuning that now must be imposed on the MSSM due to the 
LEP II constraints on the Higgs mass. 

The structure of the NMSSM can be used to motivate an extension of the GMSB messenger 
sector to include such an extra Higgs singlet field, N jl.3j. The superpotential for this sector is 
then given by: 

W = + £ N NQ$ - ^-S 2 N - ^N 2 S + \ s SH u H d + \ N NH u H d - ^N 3 . (1) 

Here S is a Standard Model gauge singlet spurion which parameterizes the supersymmetry 
breaking sector, with expectation values for the scalar and supersymmetry-breaking F compo- 
nents, (S) eM^O and (F$) = F ^ 0. The messenger fields $ and $ are in complete 5 and 5 
representations of SU(5)dSU(3)xSU(2)xU(1), so contain SU(3) triplet states with the gauge 
quantum numbers of the right handed down quarks and SU(2)-doublet states with the gauge 
quantum numbers of the left handed lepton doublet. The Higgs doublets H u and Hd couple to 
up and down type fermions, respectively. The parameters £,s,£,n,Vn, and k are expected to 
be 0(1) while rjs and Ag have to to be considerably smaller (~10 -3 ) to give a consistent phe- 
nomenology for the fi and parameters ^3| . At the tree level, the fermion components of the 



2 In some GMSB models there are additional stable SM-neutral states arising from the supersymmetry breaking 
sector, while in others the messenger and supersymmetry-breaking sectors are integrated into one sector 



leading to changes in the phenomenology |9tllCf . We will not discuss these non- minimal models in this letter. 
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$ + $ fields have, as a result of (S), a common Dirac mass, M, while because of susy-breaking 
the scalar components split into two groups with mass-squared M 2 ± F. At the one-loop level, 
the masses of the lighter group of scalars are further split, and the lightest of these scalars is the 
dark matter candidate. A computation shows that the color-charged scalars are significantly 
raised in mass, while the splitting between the lightest electrically charged, + , and neutral, 



messenger is 
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A physically acceptable model of messenger dark matter requires that this splitting be positive, 
so that the neutral scalar is the lightest state, and large enough such that charged state 4> + 
decays via + — > (fi°e + v to avoid the presence of exotic quasi-stable charged particles. 3 This 
requires — m^o > 0.6 MeV, which for M = 50 TeV imposes F/M 2 > 0.93, leading to 
m^o < 13 TeV, while for M = 150 TeV, F/M 2 > 0.96 is necessary, leading to m^o < 30 TeV. 
Thus we see that the physically relevant range of candidate dark matter messenger masses is 
-13-30 TeV. 

The potential that results from the minimizing of the superpotential Eqn. includes the 
terms: 

V = {4£ S £ N - 2£ nVn )SNM + 2 VN kSN 3 - 2 VN X N NSH u H d + ... (3) 

These three terms represent the vertices A^-l>-$, N-N-N and N-H u -Hd, respectively, each with 
a coupling enhanced by the vacuum expectation value of S, (S) ~ 100 TeV. It is because of 
this large vacuum expectation value that other terms in the potential which do not contain a 
factor of S can be safely neglected in our calculation. Furthermore, we can neglect terms that 
include the small parameters r/s and \s- 

These couplings lead to the following diagrams for messenger annihilation: 
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The thermally averaged annihilation cross sections corresponding to these diagrams are given 



by: 
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3 In Ref.|Q| it was argued that the decay had to occur before big-bang-nucleosynthesis, imposing the stringent 
constraint — m^o > 5 MeV. This is too restrictive, however, as for messenger masses in the 20 TeV 
range under consideration the total energy deposited from such decays is well below the limits set by BBN 
photodisociation. 
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and 
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where T FO is the temperature at which m^o "freezes out" of thermal equalibrium, given by 
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where Mp\ is the Planck mass, c is an 0(1) constant and g* is the number of degrees of freedom 
below the freeze out temperature (~ 232 for the particle content of the NMSSM). When solved 
by iteration, values of approximately Tpo — are found. With this in mind, we can 

see by inspecting Eqns. 0] and that the temperature independent portion of the annihilation 
cross section (the first term in each expression) dominates at the freeze-out temperature. After 
freeze-out, the relic density of <p° which remains today is given by 



Ofto/i 2 ~ 5.6 x 10 
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Thus we conclude that to obtain the cosmologically measured density of dark matter (VLp,Mh 2 ~ 
0.1), we must require (av) ~ 1.3 x 10~ 26 cm 3 /s (1.1 x 10~ 9 GeV~ 2 ). Assuming £,s,£,n,Vn, 
and k are order one constants, as they are expected to be, and (S) ~ 100 TeV, we can see 
from Eqns. 0] and that the required cross section is consistent with m^o ~ 10 TeV. In a 
parameter space scan of this model, the authors of Ref. |l3j have found masses in the range 
of approximately 7 to 24 TeV consistent with fl^oh 2 ~ 0.1. This range is somewhat artificial, 
however, depending on the range over which parameters were allowed to vary. 



III. MULTI-TEV GAMMA RAYS FROM MESSENGER DARK MATTER 

In a previous study of the HESS galactic center source jH , it was found that the a WIMP with 
a mass in the range of 12 to 29 TeV could reproduce the observed spectrum, with the best fit 
occurring for 19 TeV. This conclusion depends, however, on the assumptions made regarding 
the annihilation products of the WIMP. The results of Ref. || are therefore not completely 
general. In the case we are studying here, the dark matter candidate does not annihilate 
directly to quark or gauge boson pairs, but rather to pairs of Higgs bosons. Although the 
dominant decays of these Higgs bosons depend on the parameters chosen in this sector, in most 
models, neutral Higgses decay to either heavy quarks or gauge bosons. Therefore, the gamma- 
ray spectrum produced in the annihilations of 0°'s resembles the spectrum found in Ref. 0, 
but for a WIMP with half the mass. 

0° annihilations can also produce pairs of charged Higgs bosons. Although the gamma- 
rays produced in the decays of charged Higgses contribute minimally to the overall spectrum, 
annihilation diagrams with H + H~^ final states can contribute substantially at high energies. 
Such final states are guaranteed to be present, with the fraction of charged final states that 
include a radiated photon given by 

dN 1 d(aH+H--yv) I 'dx a (x 2 — 2x + 2) 
dx vh+h-v tt x 



m 4P 
m 2 TT+ 



x 



(8) 
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FIG. 1: The gamma-ray spectrum predicted from the annihilation of messenger dark matter particles 
(0°'s) compared to the HESS data 0|. The spectra shown represent the spectra for annihilations into 
100% neutral Higgses (solid lines) and 50% neutral and 50% charged Higgses (dashed lines). The left 
and right frames are for 20 and 25 TeV dark matter states, respectively. The results for normalized 
to best fit the HESS data. See the text for more information. 

where x = E^/m^o. As we are considering m^o ~ 20 TeV and of a few hundred GeV, 
this contribution can be significant 4 . Of course this depends on the fraction of 0° annihila- 
tions that go to charged Higgs pairs. From Eqns. H] and we can see that annihilations to 
either NN or H u Hd can dominate, for appropriate choices of parameters. For annihilations to 
H u Hdi approximately half of which produce H + H~ final states, we can expect an important 
contribution from final state photons. Annihilations to NN do not produce charged Higgses, 
as they mix only into neutral states. In figure Q we show the spectrum of gamma-rays pre- 
dicted for annihilations of 0°'s without a charged Higgs contribution (solid) and with 50% of 
annihilations producing a charged Higgs pair (dashed). We have used m h = 120 GeV and 
ttla ~ rrtH ~ vn,n± w 300 GeV, although our results depend very weakly on these choices. The 
spectra were generated using PYTHIA as implemented in the DarkSusy package 

The gamma-ray spectra predicted from messenger dark matter annihilations can provide 
good fits to the HESS data. For the spectra shown in figure [H with m^o = 25 TeV, we find 
fits with x 2 — 1-44 (solid line) and 1.20 (dashed line) per degree of freedom. For a somewhat 
lighter WIMP with m^o = 20 TeV, the quality drops to x 2 — 1-75 (solid line) and 1.36 (dashed 
line) per degree of freedom. For a 30 TeV WIMP, the x 2 's improve to 1.35 and 1.16 per degree 
of freedom with and without a charged Higgs component, respectively. It is clear that the 
presence of final state photons from charged Higgs diagrams improves the fit to the HESS data 
considerably. 

In all cases, the overall normalization of the spectrum was selected to provide the best 
possible fit to the HESS data. 

To normalize the gamma-ray spectrum from 0° annihilations, the distribution of dark matter 
near the galactic center must be considered. Unfortunately, the central region of the galaxy 
is dominated by baryons, so observations of rotation curves reveal little about the distribution 



4 Final state photons also play an important role for annihilating Kaluza-Klein dark matter which annihilate 
primarily to charge lepton pairs. See Ref. Q 
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of dark matter in this region. Furthermore, simulations have limited resolution and have 
difficulties predicting the structure of dark matter halos below the kiloparsec scale. 

With these reservations in mind, we will discuss some possibilities for the dark matter 
distribution in the central region of our galaxy. Perhaps the most popular class of distributions 
are those found by groups using N-body simulations. These include the well known Navarro 
Frenk and White (NFW) |l7| and Moore et al. |l8| profiles. These distributions display a dark 
matter cusp in the central kiloparsecs of the galaxy with the behavior p oc 1/r 7 with 7 ~ 1 
and 1.5 for the NFW and Moore et al. profiles, respectively. These N-body simulations do 
not include the effects of baryons on the dark matter distribution, however. In particular, it 
has recently been argued that the cooling of baryons in the inner kiloparsecs of our galaxy 
will compress the dark matter distribution, resulting in considerably higher densities of dark 
matter [19]. The conclusions drawn regarding this process, called adiabatic compression, are 
still somewhat controversial, however. Thirdly, the presence of a 2.6 x 10 6 solar mass black 
hole present at the dynamical center of our galaxy has lead some to argue that the accretion 



of dark matter onto this body could generate a density spike with p oc 1/r 2,4 |20j, |21[. Others 
have challenged this conclusion, however 22l |. 



The gamma-ray flux produced in the annihilations of dark matter particles near the galactic 
center is given by 

^,E,) = {av) v J N \ — 2 [ dsp 2 (r), (9) 



dE^ Anm^o Jios 

where ip is the angle from the galactic center, (crv) v ^ is the WIMP's annihilation cross section 
in the non-relativistic limit, dN 1 /dE 1 is the differential spectrum of gamma-rays produced in 
each annihilation and p(r) is the dark matter density at a distance, r, from the galactic center. 
The integral is performed over the line-of-sight of the observation. 

We can break this expression into two components, one of which depends only on the proper- 
ties of the dark matter particle and another which depends only on the dark matter distribution. 
The second of these factors is given by 

JM>) = \ ( X [ dsp 2 (r). (10) 

VVV 8.5kpcV0.3GeV/cm 3 / Ls HKJ K J 

We can then write 



$Jip, £ 7 )^6x 10- n ^ x / ,r ,, J (AQ) Aft 

7V ^' 77 d£ 7 U.3 x 10- 26 cm 3 /s/ V m o J K ' 



where J (An) is the average of J(ip) over the solid angle Aft (centered on tp — 0). Considering 
the range of masses and annihilation modes discussed in this letter, and a solid angle consistent 
with the upper limit reported by HESS (Aft ~ 3 x 10~ 6 steradians), we estimate approximately 
3 x 10 8 to be the value of J(3 x 10 _6 sr) required to generate the observed flux. The profiles 
favored by N-body simulations predict considerably smaller values, with J(3 x 10~ 6 sr) ~ 3 x 10 4 
and 4 x 10 7 for the NFW and Moore et al. profiles, respectively. These could potentially be 
increased by a factor of 10 to 100 due to adiabatic compression, however, reaching the desired 
range. The presence of a dark matter density spike from adiabatic accretion onto the central 
supermassive black hole could also potentially generate such large values. 
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IV. CONCLUSIONS 



The origin of the recently observed multi-TeV gamma-rays from the galactic center region is 
not presently understood. Although astrophysical possibilities have been proposed, annihilating 
dark matter is another possible explanation. The most popular dark matter candidates, such 
as the lightest Neutralino in models of Supersymmetry, are far too light too account for the 
very high-energy spectrum observed. 

In this letter, we have considered the possibility that stable particles from the messenger 
sector of Gauge Mediated Supersymmetry Breaking (GMSB) models are the major component 
of the dark matter of the universe and that the multi-TeV gamma-ray spectrum observed by 
HESS is the product of their annihilations. We find that 20 to 30 TeV messenger dark matter 
matter particles can reproduce the observed spectrum. The overall flux can be produced only 
if the density of dark matter is very large in the inner several parsecs of the galaxy, perhaps as 
the result of adibatic compression of the halo due to baryon cooling or the adibatic accretion 
of dark matter onto the central supermassive black hole. 

Acknowledgments: We would like to thank Joe Silk for insightful discussions. DH is sup- 
ported by the Leverhulme Trust. 
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